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ABSTRACT 

A classification scheme for close pairs of galaxies is proposed. The scheme is 
motivated by the fact that the majority of apparent close pairs are in fact wide 
pairs in three-dimensional space. This is demonstrated by means of numerical 
simulations of random samples of binary galaxies and the scrutiny of the resulting 
projected and spatial separation distributions. 

Observational strategies for classifying close pairs according to the scheme 
are suggested. As a result, physical (i.e., bound and spatially) close pairs are 
identified. 

Subject headings: galaxies: interactions — galaxies: kinematics and dynamics - 
galaxies: structure 



1. Introduction 

The investigation of binary galaxies is jeopardized by the intrinsic lack of temporal 
observational tracking of the orbital parameter space caused by orbital periods being of the 
order of hundreds of megayears. The alternative is the statistical study of samples of pairs 
relying on reasonable frequency distributions of the relevant orbital p arame ters. This has 



been done often since the pioneering work in the field by Holmberg (119371 ). Besides the 
unknowns mentioned above there is a second major setback: binary galaxy catalogues are 
always contaminated by the so-called optical pairs, i.e., unphysical, unbound pairs seen in 
projection on the plane of the sky. Earlier works, including Holmberg's, dealt with the 
problem by imposing restrictive selection criteria in the catalogues (e.g., Page et al. 1961; 
Karachentsev 1972, 1987; Turner 1976). Consequently the projected linear separations of 
pairs in these catalogues are typically 50 kpc. If one intends to study the distribution of 
galaxy mass to larger spatial extents such samples are clearly not adequate. A large sample of 
wide pairs is important in the determination of the size of dark halos, believed to exist within 
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and around the visible parts of galaxies. Later attempts have been made in this direction 
(van Moorsel 1982; Schweizer 1987; Soares 1989; Charlton & Salpeter 1991; Chengalur et al. 
1993; Nordgren et al. 1998) with the addition of wide pairs with separations as large as 1 
Mpc, and even more. The contamination by optical pairs remained a fundamental issue in 
all of these works. 

Here I focus on a definite group of pairs, namely, close pairs, to investigate the important 
point, often neglected, that closeness in the plane of the sky is not always a guarantee of 
three-dimensional closeness. 

An example from our backyard is useful as an illustration. Let us consider the most 
trivial galaxy pair: the Milky Way- Andromeda system (MW-And). Being the dominant 
galaxies in the Local Group, these galaxies form a wide pair with 700 kpc of separation. The 
cartoon in Figure [T] shows the MW-And pair and a line of sight pointing to the outskirts 
of the Virgo Cluster, located at about 15 Mpc. An observer at that location — assuming 
the Andromeda galaxy has the same linear dimension as the Milky Way — would notice a 
difference of less than 5% in their angular size. He could well classify the pair as a close one. 
With a velocity difference of 119 km s _1 at the line of sight considered, the MW-And pair 
would be a strong candidate to be in a list of close pairs. The reality is of another nature: 
MW-And is a wide pair. Its apogalacticon separation might be estimated at about 2 Mpc, if 
the system is on a highly-eccentric orbit and the galaxies have approximately 10 12 M & . Most 
plausible, the system has recently reached its apocentric configuration in a slightly-eccentric 
orbit, implying a significant transverse velocity at present. The pr oper motion of M31 has 



indeed been constrained to about 100 km s 1 by Loeb et al. (120051 ). 
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Fig. 1. — MW-And system and a line of sight towards a hypothetical distant observer. 
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Recent papers aimed at investigating a number of properties of close pairs, such as star 
formation (Barton et al. 2000; Woods et al. 2006) and merger rates, do not give due relevance 
to the fact that projected physical separation (R p ) cannot suffice as a unique constraint for 
closeness. Patton et al. (2000, 2002) define a close pair as having 5h _1 < R p < 20h _1 kpc 
and an upper limit in the pair radial- velocity difference of 500 km s _1 . The velocity criterion 
represents a somewhat loose general constraint on physical association. They are followed by 
others in essentially the same procedure (e.g., Nikolic et al. 2004). I argue here that this is 
just the first step in the selection procedure. It has to be complemented by an investigation 
of the tidal effects on the individual galaxies in the pair in order to have a clear definition 
of closeness. 

O ther s tudies address the presence of Markarian activity in galaxy pairs. Byrd & Val- 



tonen (120011 ) find that Markarian activity is an indication of tidal interaction. Their results 
should be confronted with the closeness of the sample pairs identified by an activity-blind 
procedure, as suggested here. 

In section 2 I show, using Monte Carlo simulations, that a small projected separation is 
not a secure and fair indication of a small spatial separation. This is a consequence of the 
fact that two-body orbits have most of their phase-space time at apocentric configurations 
(cf. Kepler's law). The scene is thus set for the presentation of the classification scheme 
in section 3. The relevant observational data for sorting apparently close pairs throughout 
classes are suggested in section 4. 

Section 5 summarizes the main conclusions. 



2. Spatial and projected separations in close pairs 

The most simple way of disclosing the three-dimensionality of projected separations in 
binary galaxies is by means of Monte Carlo-type simulations of bound pairs. These are 
generally an useful tool when one wants to analyze an observed binary sample in order to 
extract information such as galaxy mass and angular momentum and orbital period and 
eccentricity. Such a statistical method has been used intensively in the past (e.g., Turner 
1976; van Moorsel 1982; Schweizer 1987; Oosterloo 1988; Soares 1990, 1996). 

To simulate a sample of binary galaxies a number of input orbital parameters have 
to be specified, namely, the distribution of spatial separations and of orbital eccentricities. 
Furthermore, a mass model for galaxies has to be defined, as do projection rules to transform 
orbital velocity and radial separation into observed quantities such as line-of-sight velocity 
difference and apparent projected separation. 
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2.1. Simulation procedure 

The simulation procedure is established according to the objective one has in mind. 
Here the main goal is to investigate how close pairs, defined by their apparent separations, 
are distributed in three dimensions or, in other words, are distributed in the apparent versus 
spatial separation plane. The simulations are case problems that span a wide range of 
possibilities in the whole phase-space of binary orbits. To reach the goal I generated samples 
of 2000 artificial pairs with the following specifications. 

(i) Mass model: two equal-mass galaxies in Keplerian orbits. Galaxies move under their 

mutual point-mass Newtonian gravitational potential. I further add the restriction 
that galaxies never get closer than the sum of their (approximate) half-mass radii. I 
adopt a conservative value of 10 kpc for the sum. In practice this means that galaxies 
getting closer than this limiting separation rapidly merge and are therefore excluded 
from the simulated sample. Self-consistent iV-body simulations of binary-galaxy or- 
bit secular evolution have conclusively shown that such systems indeed merge on a 
time-scale much smaller than a Hubble time (e.g., Bartlett & Charlton 1995; Chan 
& Junqueira 2001). Evidently, galaxies are not point masses and may possess well- 
extended (dark) mass distributions. As long as the binary system is isolated from a 
third-body gravitational influence, and furthermore the mutual gravitational interac- 
tion is only radially dependent, the point-mass assumption is a reasonable description 
of the general consequences of orbital angular momentum conservation, which relies on 
Kepler's areal law, i.e., on the fact that galaxies spend most of their time at apocenter. 

(ii) Spatial separation: samples with either a fixed apocentric separation R apo of either 
200 or 600 kpc. These two values span fairly the average (apocentric) separations of 
binary galaxies as suggested by several works (Soares 1989; Schweizer 1987; Charlton 
& Salpeter 1991; Chengalur 1994; Nordgren 1997). A more rigorous simulation should 
consider a distribution of binary spatial separations. Earlier work (Gott & Turner 1979, 
but see Soares 1989, p. 17) suggested that binary galaxies follow, down to small scales, 
an extension of the classical galaxy two-point correlation function (Peebles 1980), but 
more recent evidence points to a multi-variable correlation function, depending on 
galaxy type, luminosity and local galaxy number density. An R apo value smaller than 
200 kpc is not considered here because, as mentioned above, such pairs merge quickly. 

(iii) Eccentricity distribution: three distributions, spanning the entire orbital phase-space. 
They are /(e) = 5(0.9), /(e) = 2e, and /(e) = 5(0). If binary galaxies were formed 
by early capture, suggested already by Holmberg (1937; see also Schweizer 1987), the 
triangular distribution /(e) = 2e would be suitable since it allows for some amount of 
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transverse orbital motion. The extreme cases of highly eccentric orbits (e = 0.9) and 
circular orbits bracket the triangular distribution and are considered for completeness 
and reveal fruitful consequences (see below). 

(iv) Projection on the sky: the normal to the orbital plane is distributed at random in 
space, that is, the orbital inclination i has a distribution f(i) oc sin(z). 

A pair orbit is simulated with prescriptions 1-4; velocities and radial separations are 
drawn at a random time instant within the orbit, and, from prescription 4, the projected 
separation onto the plane of the sky and the corresponding line-of-sight pair velocity dif- 
ference are calculated. The random choice of orbital phase implies, from Kepler's law, that 
pairs at apocenter — large separations — are naturally favored. For each pair of R apo and 
/(e) the process is repeated until a list of 2000 pairs is simulated, and therefrom I investigate 
the correlation between projected and spatial separation. 



2.2. Projected versus spatial separation 

A putative close pair is defined as a pair of galaxies having at most the fiducial projected 
separation R p ^ max = 50 kpc, which is typical in the studies mentioned above. I then ask: 
How many pairs with R p < R P) max have spatial separations R > R P) max ? 

Figure [2] shows N — R p histograms for the two cases considered here, R apo = 200 and 600 
kpc, and for the three distributions of eccentricities, /(e) = 6(0.9), /(e) = 2e and /(e) = 5(0). 
Shaded areas represent close pairs as defined above. Sample median values of R p are shown 
in the figure and in Table [H As expected, the least median Rp occurs for /(e) = 5(0.9) and 
for the triangular distribution /(e) = 2e, being essentially the same for both distributions. 
It amounts to 60% of R ap0 i meaning about 120 and 360 kpc respectively. In the simulations 
of pairs in circular orbits, the median R p is exactly 87% of R apo - It can be noted from Fig. 



Table 1: Number of simulated close pairs out of 2000 



Rapo 


/(e) = 5(0.9) 


/(e) = 2e 


/(e) = 5(0) 




Median Rp 


Median R p 


Median R p 


200 kpc 


284 


228 


51 




121 kpc 


125 kpc 


173 kpc 


600 kpc 


67 


30 


4 




360 kpc 


360 kpc 


520 kpc 
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[2] that the i? p distributions are the same irrespective of R apo for all eccentricities. Such an 
overall behavior is expected due to self-similarity, since the dynamical problem is the same 

- the sole difference on the orbit length scale. There is a slight fluctuation, however, in 
the triangular eccentricity distribution case (Fig. [2J middle panels), which is due to the 
additional degree of freedom introduced by the random choice of the pair eccentricity. But 
even here R p distributions are the same, with a noticeable decrease of pairs at small projected 
separations as compared to the high-eccentricity simulations (Fig. [2J top panels). 

The important difference between the distributions with R apo = 200 and 600 kpc is the 
percentage of close pairs (R p < 50 kpc). According to Tabled! 14% of the total population 
of pairs with /(e) = S(0.9) and R apo = 200 kpc are close pairs, a figure that decreases to 3% 
if R apo = 600 kpc. This is a ready consequence of the much wider phase space available for 
pairs with the larger R apo . 

As seen above, all simulations have large median R p , even greater than 2 x R Pi ma X , which 
shows that at least half of simulated pairs have spatial separations much greater than R p , m ax- 
This is further confirmed by the data shown in Table [H which represent the most important 
result, as far as the classification of close pairs is concerned: more than half of the simulated 
pairs with f? p < 50 kpc have three-dimensional separations greater than 50 kpc. 

If real binaries are predominantly in circular orbits, 100% of pairs have — from input 

- spatial separation of R > 50 kpc, and close pairs are thus very few, less than 3% of the 
simulated population. The firm conclusion from these results is that apparent separations are 
not bona-fide indicators of spatial closeness. Close pairs with R apo < 200 kpc, not included 
in the simulations, would not substantially alter the conclusion. 

In the next section, a classification scheme is put forward whose main purpose is to 
discriminate close pairs against their real physical dynamical character. The classification 
is a pre-requisite for any further astrophysical investigation of the properties of galaxies in 
close pairs. 

Table 2: Percentage of close pairs with R > 50 kpc 



Rapo 


/(e) = 6(0.9) 


/(e) = 2e 


/(e) = 6(0) 


200 kpc 


52.5% 


69.3% 


100% 


600 kpc 


53.7% 


73.3% 


100% 
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Fig. 2. — Distribution of projected separations R p for simulated pairs, with different ec- 
centricity distributions /(e) and apocentric separation R apo - Left: R apo = 200 kpe. Right: 
R ap o = 600 kpe. Dotted vertical line in all panels marks each sample's median R p . Shaded 
areas show the close-pair population. 
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3. The classification scheme 

The ground property in the classification scheme is the strength of tidal interaction 
between pair galaxies. Table [3] shows the four suggested classes of close pairs. The main goal 
is clearly discriminating between pairs that are interacting through tidal forces and those 
whose galaxies do not have their internal structure substantially affected by the mutual 
differential gravitational interaction. The strength of tidal activity between two extended 
mass distributions, with centers of mass separated by R is — grossly speaking — inversely 
proportional to R 3 , thus justifying the consideration of tidal activity as the crucial property 
for a classification scheme. As seen in Table [31 types CP I — mergers — and CP IV - 
unphysical pairs — are extreme types. Type I pairs are evolved type II pairs, and type IV 
pairs are misclassified close pairs, i.e., optical pairs. 

Any study of the effects of closeness on the properties of galaxy pairs must be preceded 
by an observational selection procedure in order to classify the sample under study into the 
four classes. In the next section suggestions for the required observational programs are 
given. 

4. Observational strategy for classifying close pairs 

4.1. CP I and CP IV 

The CP I types are easily sorted by looking to their optical disturbed morphology. 
Individual galaxies are hardly seen as separate entities. The CP IV types, on the other 
hand, are selected by their discordant redshifts. A very conservative lower limit to the line- 
of-sight galaxy velocity difference in the pair of 300 km s _1 may be adopted as an initial 
guess to distinguish optical pairs. This limit corresponds to the relative velocity of two 
1O 12 M (visible plus dark) galaxies on a circular orbit whose centers of mass are separated 
by 100 kpc. A detailed study of redshift asymmetries in the pair sample under investigation 

Table 3: Classification of close pairs 



Class 


Pair 


Interaction 


CP I 


merger 


strongest 


CP II 


tide-loud 


strong 


CP III 


tide-quiet 


weak 


CP IV 


optical 


weakest 
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might be u sed in order to evaluate the contamination by optical companions. Valtonen & 



Byrd (119861 ) find that physical pairs exhibit an equal number of positive or negative redshifts 



relative to the primary while that is not true for non physical optical pairs. 

The most difficult types to classify are CP II and CP III, which require a specially 
suitable observational strategy that takes into account the kinematical consequences of tidal 
forces. 



4.2. Tidal kinematics 

A straight application of the impulse approximation (Binney Sz Tremaine 1987) for tidal 
effects in high-speed encounters ("fly-bys") of extended objects leads to symmetrical pertur- 
bations. That is, two test particles, at symmetric positions with respect to the perturbed 
object center, acquire velocity increments of the same magnitude and opposite direction 
(Binney & Tremaine, eq. [7-54], p. 438). This will result in symmetric velocity profiles. A 
perturbed disk galaxy would thus exhibit no signature of tidal forces in its rotational profile. 

Although the impulse approximation has been successfully shown to be consistent with 
numerical experiments for mass and energy exchanges between interacting spherical galaxies 
(Aguilar & White 1985), it does not hold here. Bound pairs are mid way between the fly-by 
encounters in which the impulse approximation does apply and pairs that a re quickly evolving 



through a merger phase. Concerning the latter case, Mihos et al. (119981 ) made an iV-body 
merger model of NGC 7252, a galaxy with conspicuous tidal tails. The initial pre-merger 
galaxies have symmetrical rotation curves that evolve along the merger process to a final 
asymmetrical velocity profile (see their Fig. 6). Hence, one expects that at perigalacticon 
bound pairs will show slightly — but noticeably — distorted velocity profiles before they 
ultimately enter a merging phase. 



Furthermore, Barton et al. (119991 ) have shown, using iV-body simulations of galaxy 
encounters, that tidal interactions induce noticeable distortions on the rotation curves of 
spiral galaxy models (see their Fig. 1). 

The asymmetry of the velocity profile of both galaxies at closest approach will be used 
here as an indication of strong — loud — tidal forces, therefore implying spatial proximity. 
This is the rationale underlying the discussion of the kinematics of CP II and CP III classes 
below. 
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4.3. CP II and CP III 

The most confident criteria for CP II and CP III classification are: (1) broadband pho- 
tometry, (2) kinematics from optical single-slit spectroscopy, and (3) global H I spectral line 
low-resolution velocity profiles. These are intentionally designed to be low cost procedures 
for the classification. Strong tidal forces affect the outcome of all thr ee so rts of observations. 



Ideally, all three should be used, as, for example, in Marziani et al. (119991 ). They studied the 
close pair UGC 3995A+B (CPG 140, in Karachentsev 1972, 1987), whose POSS-II 5-band 
image is shown in Fig. [31 

Galaxy centers are separated by approximately 30", and the sm aller component B is seen 
front of the spiral disk of component A. Marziani et al. Jl999h proceeded to broadband 



in 



photometry decomposition and found unperturbed individual solutions for both galaxies. 
Optical spectroscopy reveals symmetrical kinematics relative to the galaxy centers. That is, 
the rotation-curve amplitudes are the same on the receding and approaching sides of both 
galaxies. The H I global profile is symmetrical as well. The optical velocity difference is 20 
km s _1 (inter nal error of 18 km s _1 ), and the H I velocity difference is 1 km s , from Sulentic 



& Arp (119831 ) . These observations suggest that the pair is a CP III, a "tide-quiet" close pair, 
seen at or close to apocentric orbital phase, near turnover. It is worthwhile mentioning that 
the orbital interpretation above diverges from Marziani et al. They claim that the pair is 
seen either before or after their closest approach. This seems unlikely, however because in 
this case tidal perturbations would be evident in the observations. The classification scheme 
relies on the assumption that tidal activity impresses asymmetrical signatures on morphology 
and kinematics. 

A counterexample to UGC 3995, namely, a "tide-loud" close pair, is SBG 249 (Soares 
et al. 1995). Long-slit spectroscopy has been performed and analyzed by Carvalho & 
Soares (2007, in preparation). The galaxies are NGC 1738 (leso-lv 552049011 and NGC 1739 
(leso-lv 5 520500). They have partially overlapped disks. Two spectra were obtained with the 
Double Spectrograph instrument mounted at the Palomar 5 m telescope, on 1998 February 
27-28. Figure H] shows an isophote map of the pair, with the 128" long slits superposed. The 
derived rotation curves are shown in Fig. [5] (both figures from Carvalho & Soares 2007, in 
preparation). 

The asymmetry is apparent in NGC 1738 and hinted at in NGC 1739. The southwestern 
side of NGC 1738 rotates almost twice as much as the northeastern side; that is, its closest 
side to NGC 1739 rotates at a greater speed than the farthest side. The northwestern side of 
NGC 1739 flattens out at about 100 km s _1 , while the southeastern side seems to be rising 
up at the last measured point. It is clear that the rotation profiles are strongly affected by 
their mutual tidal influence. On the other hand, the spectrum of the overlapped portion of 
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Fig. 3. — 5-band image from the POSS-II Digital Sky Survey of the close pair UGC 
3995 A+B. The smal ler ga laxy B is seen in front of the disk of the larger component. See 
also Marziani et al. (]1999l ). North is up, and east is left. 
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Fig. 4. — R-band isophote map of SBG 249 from Carvalho & Soares (2007, in preparation). 
NGC 1738 (ESO-LV 552490) is the background galaxy. Slits are 128" long. North is up, and 
east is left. 
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Fig. 5.— Rotation curve of NGC 1738 (left) and NGC 1739 (Carvalho & Soares 2007, in 
preparation). The southwestern side (SW) of the NGC 1738 disk and the northwestern side 
(NW) of NGC 1739 are partly overlapped. 
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the system is clearly disentangled showing that both disks rotate quite independently and, 
as can be seen in Fig. [5j in opposite directions. This also rules out classification as a CP I, 
in spite of the disturbed morphological appearance of the pair. 

The line of sight velocity difference is 60 km s _1 (internal error of 20 km s _1 ). Since 
they are low-mass galaxies, this is consistent with SBG 249 being a CP II seen at or close to 
perigalacticon in an eccentric orbit (cf. the Appendix in Soares 1996). Additional broadband 
photometry and global H I observation are planned for the pair. 

The strength of the asymmetry in the velocity profil e might depend on the disk rotation 



sense relative to the orbit sense. Howard et al. (119931 ) performed detailed simulations of 
tidally induced structures in disk galaxies. They discuss morphological disturbances and 
find that retrograde encounters are less effective in producing disturbances. If the result 
applies also to kinematical disturbances then pairs that pass close enough to be CP II - 
tide-loud pairs — might resemble CP III types if the galaxy's rotation sense is opposite to 
the orbit sense. These cases require careful analysis of both broadband photometry and pair 
kinematics as was the case for CPG 140 discussed above. 



5. Conclusion 

Monte Carlo simulations of close pairs of galaxies show that a projected separation 
restrictive criterion does not guarantee closeness in space. Accordingly, a simple classification 
scheme of close pairs is proposed and an observational strategy is suggested for sorting 
out a sample before any kind of astrophysical investigation of their properties is made. 
The classification is based on kinematical and dynamical criteria and on the pair's overall 
morphological appearance. The required observations are thus broadband photometry and 
optical and H I 21-cm line spectroscopy. The latter is of course only suitable for pairs with 
late-type galaxies. 

For a sample of close pairs, defined initially by their closeness in the plane of sky, the 
three kinds of observations are sufficient in ascribing to each pair one of the classes suggested. 
Classes are defined by the strength of present tidal activity and optical morphology. They 
are mergers (CP I), tide-loud pairs (CP II), tide-quiet pairs (CP III), and physically unbound 
optical pairs (CP IV). 
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for a careful reading of the manuscript. The anonymous referee is gratefully acknowledged for 
useful comments and suggestions to improve the text. The Second Palomar Observatory Sky 
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